Combination of multivariate data analysis and mixing modelling to assess tracer potential of contaminants of emerging concern in aquifers.

Collected evidence has shown that contaminants of emerging concern (CECs) in conjunction with more conventional tracers (major ions, nutrients, isotopes etc.) can be used to trace pollution origin in aquatic systems. However, in highly mixed aquifer systems signals obtained from conventional tracers overlap diminishing their potential to be used as tracers. In this study, we present an approach that incorporates multivariate statistical analysis (principal component analysis (PCA) and Kohonen's Self-Organizing Map method (SOM)) and mixing modelling to identify the most suitable CECs to be employed as anthropogenic tracers. The study area is located in the Besòs River Delta (Barcelona, NE Spain) and represents the highly mixed aquifer system. A one-year monthly based monitoring campaign was performed to collect the information about the concentrations of 105 CECs as well as major and minor ions in the river and along the groundwater flow. The dimensionality of the obtained dataset was reduced to 25 CECs, based on their estimated health risk effects, for multivariate data analysis. The obtained results showed the overlap of conventional tracers' signals obtained from PCA. In case of CECs, PCA revealed differences in their distributions allowing the differentiation of the roles of natural attenuation processes, local and regional flows on their occurrence in different parts of the aquifer. This was not possible to do using solely CECs' distribution profiles. SOMs provided the lacking information about the modality of the distribution of each CECs, revealing their ability to represent factors controlling the groundwater hydrochemistry, which assist in defining their tracer potential. Based on the obtained results four identified persistent CECs, two with unimodal (lamotrigine and 5-Desamino-5-oxo-lamotrigine) and two with bimodal (carbamazepine and diazepam (higher modality was not revealed)) distributions, were selected to run a mixing model to compare their tracer performance.

Contaminants of emerging concern in the urban aquifers of Barcelona: Do they hamper the use of groundwater?

Urban aquifers are an alternative to obtain freshwater, but they are frequently polluted by contaminants of emerging concern (CECs). Therefore, there is a need to ascertain whether CECs are a water management challenge as they might limit the use of groundwater as safe drinking water even at ng L-1 concentration levels. To answer this question, it is required to evaluate human health-risk effects of measured CECs in the groundwater and to understand their behaviour at a field-scale. This study compiles data about the presence of CECs in the aquifers of Barcelona and its metropolitan area, evaluates health risk effects of measured CECs in the groundwater and presents approaches implemented to identify and quantify the coupled hydro-thermo-chemical processes that govern their fate in the subsurface. Some CECs might be harmful to humans, such as 5-methyl-1H-benzotriazole and the pharmaceuticals azithromycin valsartan, valsartan acid, lamotrigine, gabapentin, venlafaxine and lidocaine, which show very high to intermediate health risk effects. The number of harmful CECs and the level of their hazard increase from the groups of adults and 14-18 years old teens to the groups of 4-8 years old and 1-2 years old children. Thus, some CECs can limit the use of groundwater in Barcelona as potential drinking water source. Finally, knowledge gaps in understanding the integration of these processes into urban water resources management plans are identified, which will help to define groundwater potential uses and to assure the adequate protection of the human health and the environment.

Potential of low-enthalpy geothermal energy to degrade organic contaminants of emerging concern in urban groundwater.

Low-enthalpy geothermal energy (LEGE) is a carbon-free and renewable source to provide cooling and heating to infrastructures (e.g. buildings) by exchanging their temperature with that of the ground. The exchange of temperature modifies the groundwater temperature around LEGE installations, which may contribute to enhancing the capacity of aquifers to degrade organic contaminants of emerging concern (OCECs), whose presence is significantly increasing in urban aquifers. Here, we investigate the impact of LEGE on OCECs and their bioremediation potential through numerical modelling of synthetic and real-based cases. Simulation results demonstrate that: (i) LEGE facilities have the potential to noticeably modify the concentrations of OCECs; and (ii) the final impact depends on the design of the facility. This study suggests that optimized LEGE facility designs could contribute to the degradation of OCECs present in urban aquifers, thus improving groundwater quality and increasing its availability in urban areas.

An advection-reaction model of emerging contaminants in ground and surface waterbodies with known topology: a data-based approach.

A simple data-based advection-reaction (reactive transport) model applicable to both rivers and aquifers monitoring networks is proposed. It is built on (a) available monitoring data, and (b) graph-theoretical concepts, specifically making use of the Laplacian matrix to capture the network topology and the advection process. The method yields useful information regarding the dynamic spatial behavior of the variables monitored, expressed in terms of quantitative parameters like characteristic length, entropy, first-order decay constants, synchronization between sites, and the external inputs/outputs to the system. The model was tested in an unconfined shallow aquifer located in the lower Besòs River (Spain), in which 37 pharmaceutical compounds were monitored at 7 sites, alongside two campaigns (February and May 2021). Characteristic lengths were, on average, of the same order (24.5 m) as the mean distance between consecutive monitoring sites (33.6 m), thus reflecting an adequate monitoring network design. From an estimated mean advection velocity (0.24 m·h-1), first-order decay constants were calculated for each compound and campaign, with mean values of 0.025 h-1 (February) and 0.005 h-1 (May). Whereas entropy was generally slightly larger values in February than in May (mean values of 1.02 and 0.9 entropy units respectively), synchronization showed the opposite trend (mean values of 62.4% and 68.8% respectively). The input/output profiles were generally site-dependent, regardless of the compound, and campaign considered. • A new advection-reaction modeling approach directly based on experimental data obtained from monitoring campaigns together with the network topology is proposed. • The method yields new quantitative information regarding the dynamic behavior of the variables monitored, useful for both research and management purposes.

Ecological Sustainability Assessment of Water Distribution for the Maintenance of Ecosystems, their Services and Biodiversity.

Water provision and distribution are subject to conflicts between users worldwide, with agriculture as a major driver of discords. Water sensitive ecosystems and their services are often impaired by man-made water shortage. Nevertheless, they are not sufficiently included in sustainability or risk assessments and neglected when it comes to distribution of available water resources. The herein presented contribution to the Sustainable Development Goals Clean Water and Sanitation (SDG 6) and Life on Land (SDG 15) is the Ecological Sustainability Assessment of Water distribution (ESAW-tool). The ESAW-tool introduces a watershed sustainability assessment that evaluates the sustainability of the water supply-demand ratio on basin level, where domestic water use and the water requirements of ecosystems are considered as most important water users. An ecological risk assessment estimates potential impacts of agricultural depletion of renewable water resources on (ground)water-dependent ecosystems. The ESAW-tool works in standard GIS applications and is applicable in basins worldwide with a set of broadly available input data. The ESAW-tool is tested in the Danube river basin through combination of high-resolution hydro-agroecological model data (hydrological land surface process model PROMET and groundwater model OpenGeoSys) and further freely available data (water use, biodiversity and wetlands maps). Based on the results, measures for more sustainable water management can be deduced, such as increase of rainfed agriculture near vulnerable ecosystems or change of certain crops. The tool can support decision making of authorities from local to national level as well as private enterprises who want to improve the sustainability of their supply chains.

Cooling-induced reactivation of distant faults during long-term geothermal energy production in hot sedimentary aquifers.

Deep geothermal energy (DGE) represents an opportunity for a sustainable and carbon-free energy supply. One of the main concerns of DGE is induced seismicity that may produce damaging earthquakes, challenging its widespread exploitation. It is widely believed that the seismicity risk can be controlled by using doublet systems circulating water to minimize the injection-induced pressure changes. However, cold water reinjection may also give rise to thermal stresses within and beyond the cooled region, whose potential impacts on fault reactivation are less well understood. Here, we investigate by coupled thermo-hydro-mechanical modeling the processes that may lead to fault reactivation in a hot sedimentary aquifer (HSA) in which water is circulated through a doublet. We show that thermal stresses are transmitted much ahead of the cooled region and are likely to destabilize faults located far away from the doublet. Meanwhile, the fault permeability mainly controls the fault reactivation timing, which entails the importance of employing appropriate characterization methods. This investigation is crucial for understanding the mechanisms controlling induced seismicity associated with DGE in a HSA and allows the success of future DGE projects.

Fate and risk assessment of sulfonamides and metabolites in urban groundwater.

Antibiotics, such as sulfonamides (SAs), have recently raised concern as wastewater treatment plants (WWTPs) partly remove them, and thus, SAs continuously enter the aquifers. In this context, the aims of this work are to (1) investigate the temporal evolution of SAs and metabolites in an urban aquifer recharged by a polluted river; (2) identify the potential geochemical processes that might affect SAs in the river-groundwater interface and (3) evaluate the ecological and human health risk assessment of SAs. To this end, 14 SAs and 4 metabolites were analyzed in river and urban groundwater from the metropolitan area of Barcelona (NE, Spain) in three different sampling campaigns. These substances had a distinct behavior when river water, which is the main recharge source, infiltrates the aquifer. Mixing of the river water recharge into the aquifer drives several redox reactions such as aerobic respiration and denitrification. This reducing character of the aquifer seemed to favor the natural attenuation of some SAs as sulfamethoxazole, sulfapyridine, and sulfamethizole. However, most of the SAs detected were not likely to undergo degradation and adsorption because their concentrations were constant along groundwater flow path. In fact, the intensity of SAs adsorption is low as the retardation factors are close to 1 at average groundwater pH of 7.2 for most SAs. Finally, risk quotients (RQs) are used to evaluate the ecological and human health risks posed by single and mixture of SAs in river water and groundwater, respectively. Life-stage RQs of the SAs detected in groundwater for the 8 age intervals were low, indicating that SAs and their mixture do not pose any risk to human beings. Concerning the environmental risk assessment, SAs do not pose any risk for algae, fish and crustaceans as the RQs evaluated are further lower than 0.1.

Occurrence of pathogens in the river-groundwater interface in a losing river stretch (Besòs River Delta, Spain).

The aim of this study is to investigate the occurrence of faecal indicator and microbial pathogens (bacteria and virus) in the shallow urban aquifer of the Besòs River Delta (NE Spain). To this end, human adenovirus (HAdV) and Norovirus of genogroups I and II (NoV GI and NoV GII) as well as the faecal indicator bacteria (FIB) Escherichia coli (EC) and faecal enterococci (FE) were monitored in groundwater and in the River Besòs in December 2013 and in July 2104. None of the targeted pathogens were detected in groundwater in December 2013 but contamination of human origin was observed in approximately 50% of the points sampled in July 2014 reaching concentrations up to 99 GC/100 mL for HAdV. Generally, microbial concentrations in river water were higher than those detected in groundwater. This observation indicates that pathogens are naturally attenuated when river water infiltrates and flows through the aquifer, however HAdV were detected at a sampling point located at 380 m from the river in the absence of FIB. The presence of human viral contamination may represent a risk for the use of groundwater as a drinking water source. Further research is needed to understand the dynamics of pathogens in river-groundwater interface over long time periods and a wide range of flow conditions (wet and dry periods) since the urban groundwater of this aquifer might be a valuable drinking water resource in Barcelona especially during drought periods. The methodology followed in this research can be applied to other urban aquifers with similar purposes since the scarcity and contamination of freshwater resources are worldwide issues.

Parametric assessment of hydrochemical changes associated to underground pumped hydropower storage.

Underground pumped hydropower storage (UPHS) using abandoned mines is an alternative to store and produce electricity in flat regions. Excess of electricity is stored in form of potential energy by pumping mine water to a surface reservoir. When the demand of electricity increases, water is discharged into the mine (i.e., underground reservoir) through turbines producing electricity. During the complete operational process of UPHS plants, hydrochemical characteristics of water evolve continuously to be in equilibrium successively with the atmosphere (in the surface reservoir) and the surrounding porous medium (in the underground reservoir). It may lead to precipitation and/or dissolution of minerals and their associated consequences, such as pH variations. Induced hydrochemical changes may have an impact on the environment and/or the efficiency (e.g., corrosions and incrustations affect facilities) of UPHS plants. The nature of the hydrochemical changes is controlled by the specific chemical characteristics of the surrounding porous medium. However, the magnitude of the changes also depends on other variables, such as hydraulic parameters. The role of these parameters is established to define screening criteria and improve the selection procedure of abandoned mines for constructing UPHS plants. This work evaluates the role of the main hydrogeological factors for three different chemical composition of the porous medium. Results are obtained by means of numerical reactive transport modeling. Potential impacts on the environment (mainly on groundwater and surface water bodies) and on the efficiency of the UPHS plants vary considerably from a hydraulic parameter to another showing the need for a detailed characterization before choosing locations of future UPHS plants.

Occurrence of greenhouse gases in the aquifers of the Walloon Region (Belgium).

This work aims to (1) identify the most conductive conditions for the generation of greenhouses gases (GHGs) in groundwater (e.g., hydrogeological contexts and geochemical processes) and (2) evaluate the indirect emissions of GHGs from groundwater at a regional scale in Wallonia (Belgium). To this end, nitrous oxide (N2O), methane (CH4) and carbon dioxide (CO2) concentrations and the stable isotopes of nitrate (NO3-) and sulphate were monitored in 12 aquifers of the Walloon Region (Belgium). The concentrations of GHGs range from 0.05µg/L to 1631.2µg/L for N2O, 0µg/L to 17.1µg/L for CH4, and 1769 to 100,514ppm for the partial pressure of CO2 (pCO2). The highest average concentrations of N2O and pCO2 are found in a chalky aquifer. The coupled use of statistical techniques and stable isotopes is a useful approach to identify the geochemical conditions that control the occurrence of GHGs in the aquifers of the Walloon Region. The accumulation of N2O is most likely due to nitrification (high concentrations of dissolved oxygen and NO3- and null concentrations of ammonium) and, to a lesser extent, initial denitrification in a few sampling locations (medium concentrations of dissolved oxygen and NO3-). The oxic character found in groundwater is not prone to the accumulation of CH4 in Walloon aquifers. Nevertheless, groundwater is oversaturated with GHGs with respect to atmospheric equilibrium (especially for N2O and pCO2); the fluxes of N2O (0.32kgN2O-NHa-1y-1) and CO2 (27kgCO2Ha-1y-1) from groundwater are much lower than the direct emissions of N2O from agricultural soils and fossil-fuel-related CO2 emissions. Thus, indirect GHG emissions from the aquifers of the Walloon Region are likely to be a minor contributor to atmospheric GHG emissions, but their quantification would help to better constrain the nitrogen and carbon budgets.

Occurrence, fate and risk assessment of personal care products in river-groundwater interface.

This work presents the occurrence and fate of selected personal care products (PCPs) in the urban river-groundwater interface. To this end, urban river and groundwater samples were collected in Sant Adrià del Besòs (NE of Spain) and a total of 16 PCPs were analyzed including benzophenone derivatives, camphor derivatives, p-aminobenzoic acid derivatives, triazoles and parabens in three different campaigns (from May 2010 to July 2014). These compounds reach the aquifer through the recharge of Besòs River that receives large amounts of effluents from waste water treatment plants. Results have shown that most of the compounds were not or barely detected (maximum concentrations around 200ng/L) in groundwater samples during the different sampling campaigns. Only two triazoles, namely benzotriazole (BZT) and methyl benzotriazol (MeBZT) were found at high concentrations in groundwater samples (maximum concentration around 2000ng/L). The fate of PCPs in the aquifer was assessed using mixing analysis considering the seasonal variability of the Besòs River. Overall, measured groundwater concentrations were significantly much lower than those estimated by the mixing of the river water. This observation suggested that most of the PCPs are naturally removed when river water infiltrates the aquifer. However, some compounds were more persistent in the aquifer. These compounds were in descending order: the triazoles BZT and MeBZT followed by the camphor derivative 4MBC and the paraben MePB. The measured concentrations allowed us to assess the environmental risk posed by the selected UV-filters and parabens in the river and groundwater samples. Hazard Quotients (HQs) for different aquatic species were calculated in order to characterize the ecotoxicity potential of the studied compounds in the river-groundwater interface. HQ values were always below 1 indicating that at the concentrations observed in the surface or aquifer water of Besòs River these compounds pose no risk to the selected aquatic organisms.

Quantifying chemical reactions by using mixing analysis.

This work is motivated by a sound understanding of the chemical processes that affect the organic pollutants in an urban aquifer. We propose an approach to quantify such processes using mixing calculations. The methodology consists of the following steps: (1) identification of the recharge sources (end-members) and selection of the species (conservative and non-conservative) to be used, (2) identification of the chemical processes and (3) evaluation of mixing ratios including the chemical processes. This methodology has been applied in the Besòs River Delta (NE Barcelona, Spain), where the River Besòs is the main aquifer recharge source. A total number of 51 groundwater samples were collected from July 2007 to May 2010 during four field campaigns. Three river end-members were necessary to explain the temporal variability of the River Besòs: one river end-member is from the wet periods (W1) and two are from dry periods (D1 and D2). This methodology has proved to be useful not only to compute the mixing ratios but also to quantify processes such as calcite and magnesite dissolution, aerobic respiration and denitrification undergone at each observation point.

Urban groundwater contamination by residues of UV filters.

The occurrence and fate of UV filters (UV F) in an urban aquifer in correlation with (1) the spatial distribution of UV F in Barcelona's groundwater, (2) the depth of the groundwater sample, (3) the physicochemical properties of the target compounds, (4) the recharge sources, and (5) the redox conditions of the Barcelona aquifers, were studied for the first time. The highest groundwater concentrations and the largest number of detected UV F were observed in an aquifer recharged by a polluted river (around 55 ng/L in SAP-4). In contrast, the urbanized areas had lower concentrations (around 20 ng/L in MPSP-1). Two pathways can be identified for UV F to enter the aquifers: (1) leakage of row sewage from the sewage network in urbanized areas and (2) wastewater treatment plant (WWTP) effluents discharged into the river. Measured concentrations of UV F were significantly much lower than those estimated from the waste water proportion in groundwater samples suggesting that UV F might undergo transformation processes in both reducing and oxidizing conditions.

Occurrence of carbamazepine and five metabolites in an urban aquifer.

This paper deals with urban groundwater contaminated with carbamazepine (CBZ) and five of its human metabolites in Barcelona. Groundwater samples were accordingly collected in the aquifers of Poble Sec and Besòs River Delta. Higher concentrations and more compounds were found in the Besòs River Delta aquifer, which is recharged by a river contaminated with treated effluent from numerous treatment plants. By contrast, the urban area of Poble Sec presented lower concentrations and fewer compounds. The results showed that CBZ could be attenuated in the Poble Sec aquifer since concentrations in groundwater were lower than those evaluated from mixing of the recharge sources. Conversely, CBZ and its human metabolites were not removed under the reducing conditions of the Besòs River Delta aquifer probably because of the short residence time in this aquifer.

Emerging organic contaminants in groundwater in Spain: a review of sources, recent occurrence and fate in a European context.

This paper reviewed the presence of emerging organic contaminants (EOCs) that have been found in the groundwater in Spain in both, rural and urban areas. The list of compounds included pesticides, pharmaceutical active compounds (PhACs), selected industrial compounds, drugs of abuse (DAs), estrogens, personal care products and life-style compounds. The main sources of pollution and possible pathways have been summarised in this review. EOCs are likely to enter to the aquifer mainly through the effluents of waste water treatment plants (WWTPs) and are present in groundwater at concentrations of ng/L to µg/L. The most studied compounds in Spanish groundwater were pesticides followed by industrial compounds and PhACs. It is important to mention that compared to other water bodies, such as rivers, groundwater is considerably less contaminated, which may be indicative of the natural attenuation capacity of the aquifers. However, some EOCs have sometimes been detected at higher concentration levels in the aquifer than in the rivers, indicating the need for further research to understand their behaviour in the aquifers. For a wide array of compounds, their maximum concentrations show values above the European groundwater quality standard for individual pesticides (0.1 µg/L). Therefore, to preserve groundwater quality against deterioration it is necessary to define environmental groundwater thresholds for the non-regulated compounds.

Drugs of abuse in urban groundwater. A case study: Barcelona.

This study is concerned with drugs of abuse (DAs) and their metabolites in urban groundwater at field scale in relation to (1) the spatial distribution of DAs in Barcelona's groundwater, (2) the depth of the groundwater samples, (3) the presence of DAs in recharge sources, and (4) the assessment of the fate of DAs in Barcelona aquifers. To this end, 37 urban groundwater samples were collected in the city of Barcelona and a total of 21 drugs were analyzed including cocainics, amphetamine-like compounds, opioids, lysergics and cannabinoids and the prescribed drugs benzodiazepines. Overall, the highest groundwater concentrations (around 200 ng/L in SAP-4) and the largest number of detected DAs were found in zones basically recharged by a river that receives large amounts of effluents from waste water treatment plants (WWTPs). In contrast, the urbanized areas yielded not only lower concentrations but also a much smaller number of drugs, which suggests a local origin. In fact, cocaine and its metabolite were dominant in more prosperous neighborhoods, whereas the cheaper MDMA was the dominant DA in poorer districts. Measured concentrations were consistently smaller than those estimated from the waste water fraction in groundwater samples, suggesting that DAs undergo removal processes in both reducing and oxidizing conditions.